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Cosmological constraints on generalized Chaplygin gas model: Markov Chain Monte 
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Dalian University of Technology, Dalian, 116024, P- R- China 

We use the Markov Chain Monte Carlo method to investigate a global constraints on the gen- 
eralized Chaplygin gas (GCG) model as the unification of dark matter and dark energy from the 
' latest observational data: the Constitution dataset of type supernovae la (SNIa), the observational 

' Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), 

^ l' and the cosmic microwave background (CMB) data. In a non-fiat universe, the constraint results 

<C ; for GCG model are, Q,h^ = 0.0235t°;[JEJ?^ (la) t^oZll (2^), ilk = 0.0035^°:°}^ (la) l^^gi (2a), 

O ■ A.^ 0.753t°:«^3^ (1'-) tHS (2a), « = OM3+_°o fol i^^) t^:i?7 (2't), and Ho = rOm+JH (la) 

(2a), which is more stringent than the previous results for constraint on GCG model parameters. 
Furthermore, according to the information criterion, it seems that the current observations much 
support ACDM model relative to the GCG model. 
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Introduction 



> , 

. Recently, mounting cosmic observations suggest that the expansion of present universe is speeding up rather than 
'NT ■ 

CO ■ slowing dovifn 1]. And they indicates that baryon matter component is about 5% for total energy density, and about 

cn ■ 

• ■ 95% energy density in universe is invisible. Considering the four-dimensional standard cosmology, this accelerated 

■ expansion for universe predict that dark energy (DE) as an exotic component with negative pressure is filled in 
' universe. And it is shown that DE takes up about two-thirds of the total energy density from cosmic observations. 

^ On the other hand, in theory many kinds of DE models _have already been constructed in order to explore the DE 

; properties. For a review on DE models, please see Refs. 

■ It is well known that the generalized Chaplygin gas (GCG) model have been widely studied for interpreting the 



accelerating universe [3||[J]. The most interesting property for this scenario is that, two unknown dark sections in 
universe-dark energy and dark matter can be unified by using an exotic equation of state. In this paper, we use 
the Markov Chain Monte Carlo (MCMC) technique to constrain the GCG model from the latest observational data: 
the Constitution dataset including 397 type la supernovae (SNIa), the observational Hubble data (OHD) fl, the 
cluster X-ray gas mass fraction [7|, the measurement results of baryon acoustic oscillation (BAO) from Sloan Digital 
Sky Survey (SDSS) and Two Degree Field Galaxy Redshift Survey (2dFGRS) QlQ, and the current cosmic microwave 
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background (CMB) data from five-year WMAP 



0. 



2. generalized Chaplygin gas model 



For GCG model, the energy density p and pressure p are related by the equation of state 



a 



PGCG - -— , (1) 

Pgcg 

where A and a are parameters in the model. By using the energy momentum conservation equation in the Friedmann- 
Robertson- Walker (FRW) cosmology, d{pa^) = —pd{a?), the energy density of GCG fluid can be given 



PGCG = POGCG 



[A, + (l-A,)(l + z)3(i+")]TT^, (2) 



where a is the scale factor, z describes the redshift, and = — It is easy to see that As denotes the current 

Pa 



equation of state for GCG fluid 



Il[[l2|. Since the GCG fluids behaves as dust at early stage and as dark energy at 
later stage, it can be considered as a scenario of the unification of dark matter and dark energy [isl. Considering 
a non-flat universe is filled with three components: the GCG component, the baryon matter component, and the 
radiation component, we have the total energy density, ptotai — Pgcg + Pb + Pr + Pk- Making use of the Friedmann 
equation, the Hubble parameter H is expressed as 

= ^ -n^-Qr- ^k)[As + (1 - ^)(i + + n^ii + zf + a-(i + z)^ + n,{i + z)'}, (3) 

where Hq — 100ft. km s^^Mpc^^ is the present Hubble constant, ilb, ilr, and ftk denote dimensionless baryon matter, 
radiation, and curvature density, respectively. 



3. Current observational data and cosmological constraint 

In this section, we introduce how the currently available data are used to constrain the model in our calculation. 

3.1. Type la super novae 

We constrain the parameters with Constitution dataset Q including 397 SNIa Q, which is obtained by adding 90 
SNIa from CfA3 sample to 307 SNIa Union sample [l^. CfA3 sample are all from the low-redshift SNIa, z < 0.08, 
and these 90 SNIa are calculated with using the same Union cuts. The addition of CFA3 sample increases the number 
of nearby SNIa and reduces the statistical uncertainties. The theoretical distance modulus p{z)th is defined as 

^ith{z) = ^\og^a[DL{z)]+ (4) 

In this expression Dl{z) is the Hubble-free luminosity distance HodL{z)/c, and 

, , ^ C{l + Z) . r— r dz' 



VW\ Jo Hiz'Y 

PLQ = 42.38- 51ogio/i. 



where sinn{y^\Qk\x) respectively denotes sin(-y/|rife|x), -y/jTIfcla::, sinh(-\/|rifc|a;) for flk < 0, = and Qk > 0. 
Additionally, the observed distance moduli ^obs{zi) of SNIa at Zi is 

tJ.obs{zi) mobs{zi) - M, (5) 

where M is their absolute magnitudes. 

For the SNIa dataset, the best fit values of the parameters Ps can be determined by a likelihood analysis, based on 
the calculation of 

X [Ps,M ) = 2^ -2 

SNIa » 

^ ^ {51ogio[DL(p^,^:i)] -mofc3(zj) + M'}'^ 

SNIa 

where ps denotes the parameters contained in the model, M' = /io + Af is a nuisance parameter which includes the 
absolute magnitude and the parameter h. The nuisance parameter M' can be marginalized over analytically [fsj l as 



/-i-oo 
exp 
-OO 

resulting to 



+ 00 



with 



. _ {5\ogj^„[DL{ps,z^)]-mobs{zt)}'^ 

SNIa ' 

5log-^Q[DL{ps,Zi)] - mobsizt) 



SNIa 



SNIa 



Relation © has a minimum at the nuisance parameter value M' = B/C, which contains information of the values of 
h and M. Therefore, one can extract the values of h and M provided the knowledge of one of them. Finally, note 
that the expression 

xlMiaiPs,B/C)=A-{ByC), (8) 



which coincides to ([71) up to a constant, is often used in the likelihood analysis jl5|, and thus in this case the results 
will not be affected by a flat M' distribution. 



3.2. Observational Hubble data 



li-In 



0, 



Jimenez et al. obtained an 



16|, and used it to constrain the EOS 



The observational Hubble data are based on differential ages of the galaxies 
independent estimate for the Hubble parameter using the method developed in 
of dark energy. The Hubble parameter depending on the differential ages as a function of redshift z can be written in 
the form of 



So, once dz/dt is known, H{z) is obtained directly [6|. By using the differential ages of passively-evolving galaxies 
from the Gemini Deep Deep Survey (GDDS) [l8| and archival data 18 ; 24 L Simon et al. obtained H{z) in the 
range of < z < 1.8 The twelve observational Hubble data from [313 are list in Table ID In addition, in 



2 





0.1 


0.17 0.27 


0.4 


0.48 


0.88 


0.9 


1.30 1.43 


1.53 


1.75 


H{z) (km s-^Mpc-i) 


74.2 


69 


83 77 


95 


97 


90 


117 


168 177 


140 


202 


Icr uncertainty 


±3.6 


±12 


±8 ±14 


±17 


±60 


±40 


±23 


±17 ±18 


±14 


±40 



IQ. 



TABLE I: The observational H{z) data [2 

[2^, the authors took the BAO scale as a standard ruler in the radial direction, obtain three more additional data: 
H{z = 0.24) = 79.69 ± 2.32, H{z = 0.34) = 83.8 ± 2.96, and H{z = 0.43) = 86.45 ± 3.27. 



The best fit values of the model parameters from observational Hubble data are determined by minimizing 



28| 



15 



(10) 



where Hth is the predicted value for the Hubble parameter. Hobs is the observed value, a(zi) is the standard deviation 
measurement uncertainty, and the summation is over the 15 observational Hubble data points at redshifts Zi. 



3.3. The X-ray gas mass fraction constraints 

According to the X-ray cluster gas mass fraction observation, the baryon mass fraction in clusters of galaxies (CBF) 
can be utilized to constrain cosmological parameters. The X-ray gas mass fraction, fgas, is defined as the ratio of the 
X-ray gas mass to the total mass of a cluster, which is a constant and independent on the rcdshift. In the framework 
of the KCDM reference cosmology, the X-ray gas mass fraction is presented as Q 



J qas \Z) 









\ Da{z) \ 



(11) 



'f''^ ' ' l + s{z) 

where A is the angular correction factor, which is caused by the change in angle for the current test model ^2500 in 
comparison with that of the reference cosmology 6*2500'*'^ • 



A = 



f0Mi!^y ^ ( h{z)Da{z) y 

I 02500 ) ~ \ [H{z)Da{z)]''^'''' ) ' ^ ^ 



here, the index 77 is the slope of the fgas{r/r25oo) data within the radius ^25007 with the best-fit average value 
ry = 0.214 ± 0.022 7]. And the proper (not comoving) angular diameter distance is given by 

r- , fir' 

It is clear that this quantity is related with di^{z) by 

dL{z) 



Da{z) 



(l + z)2- 

For GCG model, since it is considered as the unification of dark matter and dark energy, we do not have dark matter 
in this model. So, the matter density is not explicitly included in the background equation Following the Ref. 



[2^, we use an relation between and Ag'. ilm = 1 — ^ ^ ^s(l — ^ — ^k)- This parameter, fim, 
is an estimate of the "matter" component of the GCG fluid with the baryon density. And for the BAO and CMB 
constraint methods in the following, we also take this expression of f2„i. 

In equation (jlip . the parameter 7 denotes permissible departures from the assumption of hydrostatic equilibrium, 
due to non-thermal pressure support; the bias factor b(z) = feo(l + C(bz) accounts for uncertainties in the cluster 
depletion factor; s(z) = so(l + agz) accounts for uncertainties of the baryonic mass fraction in stars and a Gaussian 
prior for sq is employed, with sq ~ (0.16 ± 0.05)h^Q Q; the factor K is used to describe the combined effects of the 
residual uncertainties, such as the instrumental calibration and certain X-ray modelling issues, and a Gaussian prior 
for the 'calibration' factor is considered by if = 1.0 ± 0.1 [7]; 

Following the method in Ref. 0, and adopting the updated 42 observational fgas data in Ref . Q , the best fit 
values of the model parameters for the X-ray gas mass fraction analysis are determined by minimizing, 

N rrACDMf A _ f / ^l2 
2 Unas \^tJ Jgasy^i)\ . 

. (14) 

i J gas ^ 

3.4. Baryon acoustic oscillation 

The baryon acoustic oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, and measure the distance-redshift relation at 2: = 0.2. Additionally, baryon acoustic oscillations in the 
clustering of the SDSS luminous red galaxies measure the distance-redshift relation sd z — 0.35. The observed scale 
of the BAO calculated from these samples, as well as from the combined sample, are jointly analyzed using estimates 
of the correlated errors to constrain the form of the distance measure Dv{z) [ol. I31I1 



Dv{z) = 



{l + zrD\{z)- 



1/3 



(15) 



Hiz)_ 

The peak positions of the BAO depend on the ratio of Dv{z) to the sound horizon size at the drag epoch (where 
baryons were released from photons) Zd, which can be obtained by using a fitting formula [3^: 

lOQUO 7,2^-o.4l9 

with 

bi = 0.3l3{n,nh^)-°-^^^[l + 0.607(an/^')°-^'^], (17) 
&2 = 0.238{n„,hy-^^^. (18) 

In this paper, we use the data of rs{Z(i)/ Dv{z) extracted from the Sloan Digitial Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) S], which are listed in Table HIl where rs{z) is the comoving sound 



horizon size 



, , /■* Csdt r csda c 
rs(z) =c — — = c / dz 



c 



Q a Jq a?H H{z) 
1/(1+^) da 

(19) 



where Cg is the sound speed of the photon— baryon fluid [3. 



3 P-y{z) 3 ^-y 

and here = 2.469 x IQ-^/i^^ foj. Tca/b = 2.725K. 



2 = 3+- X -^-^i^ ==3 + - X (-^)a, (20) 



z 


rs{zd)/Dv{z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE II: The observational rs{zd)/Dv{z) data 

Using the data of BAO in Table HIl and the inverse covariance matrix V^^ in 

/ 30124.1 -17226.9 \ 
V-'^\ , (21) 

\^ -17226.9 86976.6 J 

thus, the XsAoiPs) given as 

xlAoiPs)=X'V-'X, (22) 
where X is a column vector formed from the values of theory minus the corresponding observational data, with 

^^^^ - 0.190533 



^ = I "^l":?^ I - (23) 

and X* denotes its transpose. 



3.5. Cosmic microwave background 



The CMB shift parameter R is provided by 



R{z,) = ^n„,Hl{l + z,)Da{z.,)/c, (24) 

which is related to the second distance ratio DA{Zf,)H{z^/c by a factor 1 ± z*. The redshift z* (the decoupling 
epoch of photons) is obtained using the fitting function [36| 

z, = 1048 [1 + 0.00124(f76/i2)-0-738j ^ g^(n„,hy^] , 
where the functions gi and g2 read 

51 = 0.0783(f^6/l2)-0-238(^^39 5(f^^;^2)0.763)-l^ 

52 = 0.560 (1 + 21. 1(r2b/i2)isi)"\ 

In additional, the acoustic scale is related to the first distance ratio, DA{z*)/rs{z:,) (One can define an angular scale 
of the sound horizon at decoupling epoch), and is defined as 

nDAjz^) . 
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parameters flat GCG non-flat GCG flat ACDM non-flat ACDM 





n n':)'?'?+o o''23+o.oo29 

U.UZOO_g QQJg_Q QQ20 


n n9'ic:+0 0021+0.0028 
U.UZOO_Q QQj^g_Q Q022 










flu 




A nn'<c;+0 0172+0.0226 
U.UUOO_Q oj^g2_o.o204 








U.UUUZ_0 0119-0.0127 








0.2801 


0.029- 
0.030 


f0.039 
-0.038 


n 070+0.027+0.037 
U.z/O_o 031-0.039 


As 


n 7ftn+0-029+0.034 
U. ( U"_o. 039-0. 046 


n 7c;q+0.037+0.045 
U. (00_g Q35_Q 










a 


n nqq+0. 066+0. 096 
U.UJJ_0. 071-0. 087 


n n,4'^+0102+0.134 
U.U'lO.o. 106-0. 117 










Ho 


CQ Q7+2.87+3.48 
Dy.y(_2.7g_3.08 


(U.UU_2.g2_3.67 


70.01 


+2.51- 
-2.01 


f3.19 
-2.69 


70 9/1+2.55+3.36 
'"•^^-3.34-3.90 



xLn(Xmin/do/) 519.342 (1.139) 519.371 (1.144) 520.351 (1.139) 520.302 (1.141) 



TABLE III: The data fltting results of GCG and ACDM model parameters with la and 2a confidence levels for flat and non-flat 
universe. 



where a factor of 1 + arises because Da{z) is the proper (physical) angular diameter distance, whereas rs(z,) is the 
comoving sound horizon at z*. Using the data of I a, R, z^, in 10[ and their covariance matrix of [Z^(z,), R{z^,), z-^,], we 
can calculate the likelihood L as Xc'J\/b ^ ^ 

xImb = l^d,[Cov-\d,,d,)[l^d^]% (26) 

where Ad^ — di — df°'^°' is a row vector, and di — [Ia, R,z^). 



4. Method and results 

According to the descriptions in section[3j for our calculations the total likelihood function is written as L cx e~^^/^, 
here the equals 

— XSNIa + XoHD + XcBF + XbAO + XcMB- (27) 

In our analysis, we perform a global fitting on determining the cosmological parameters using a MCMC. The MCMC 
code is listed in the publicly available CosmoMC package [s^ written in Fortran 90. In addition, the likelihood of 
fgas has been included in the modified CosmoMC. For the analysis of X-ray cluster gas mass fraction, we seek help 
from the online Fortran 90 code 0, ^^.d correct the patch for fgas after some nontrivial crosschecks. For each 

MCMC calculation on GCG model, we run 8 independent chains comprising of 50000-60000 chain elements. The 
average acceptance rate is about 35%. To get the converged results, we test the convergence of the chains by typically 
getting i? — 1 to be less than 0.03. 

In Fig. [U we show one dimensional probability distribution of each parameter and two dimensional plots for 
parameters between each other for the GCG model in flat and non-flat universe. Dotted lines are mean likelihoods 
of samples, and solid lines are marginalized probabilities for ID distribution in the figure. According to Fig. [1] the 
constraint results on the best fit values of cosmological parameters with la and 2a confidence levels are listed in 
Table IIIIl From this table, we can see that this constraint on GCG model parameter, Q^h^ — 0.0233lQQQ^g (Ict) 
-ofolo (2'^), ^ = Q-760-om9 (1^^2,046 (2c^) and a = 0.033^°;°?^ (la) tom7 (2^^) for a flat universe is more 



stringent than the result in Refs. 



40|. Furthermore, we also consider the constraint on GCG model in a non-flat 
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FIG. 1: 1-D constraints on model parameters and 2-D contours on these parameters with la, 2a confidence levels between each 
other in flat (left) and non-flat (right) GCG model. 




26 3 68 72 -0.01 0.01 0.25 0.3 66 68 70 72 

Ho "k "o 



FIG. 2: 1-D constraints on model parameters and 2-D contours on these parameters with la, 2a confidence levels between each 
other in flat (left) and non-flat (right) ACDM model. 

universe, and the constraint result are: rjfc/i^ = OmSbtomil (lo") ^0:0022 (2cr), = 0-0035t°o °olll (Icr) tomol (2'^), 
As = 0.753toml (1^) (2cr), and a = OMSto.wl i^^) -o.iiv (2^). From our results above, it is clear that 

the best fit parameter a has a very small value, and its value is a little bigger than zero. Furthermore, for the case 
of a = 0, it is included in the la confidence level around the best fit parameter. One knows for a = 0, the GCG 
model reduces to the cosmic concordance model, ACDM. Thus, it seems that the current cosmic observations indicate 
the concordance model should be favored when we apply the present observed data to constrain a more complex 
GCG model. Then in the following, we also consider the constraint on ACDM model from the above datasets, and 
the results are: n„, = 0.280t^;J^3^ (Ict) to fal {2a), Ha = 70.0ll2:oi (Icr) ^2:69 (2'^) with x™,„ = 520.351 for a flat 
universe, and = 0.0002t°;™83 (1^) +o.oi46 (2^)^ ^ 0.278t°:°27 (la) ^^011 (2^^), Ho = 70.24^2-535 (i^) +3-36 
{2a) with Xmin — 520.302 for a non-flat case. It can be seen that though ilk is not exactly equal to zero, the curvature 
density parameter has a very small value, and the best fit result shows a close universe. 

We also calculate the values of Xmin/'^"/ for GCG and ACDM model in Table ITTll where the value of dof (degree of 
freedom) equals to the number of observational data points minus the number of parameters. From Table Hill we can 
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Case model 


Free model parameters 


2 


AAIC 


ABIC 


flat ACDM 


f^m, Ho 


520.351 








non-flat ACDM 


^k, f2m, Ho 


520.302 


1.951 


6.080 


flat GCG 


Qb, As, a. Ho 


519.342 


2.991 


11.249 


non-flat GCG 


flk, ^b. As, a. Ho 


519.371 


5.020 


17.407 



TABLE IV: Information criteria results 



see that for using the large data points, the difference of Xmin/'^'^f ^'^^ ^^'^ different models is not obvious. It appears 
that for applying the quantity X^inl dof to rate the goodness of models, it is not a good (or refined) method. So, in 
the foUowing we also use the objective information criteria (IC) to estimate the quality of the models. The Akaike 
information criteria (AIC) and Bayesian information criteria (BIC) are respectively defined as 

AIC = -2\nCmax + '2K, (28) 



BIC = -2\nC,nax + K\nn, (29) 

where Umax is the highest likelihood in the model with the best fit parameters, K is the number of estimable 
parameters, n is the number of data points in the fit^. The term, — 21n£,„aa; = Xmim measures the quality of model 
fit, while the terms including K interpret model complexity. The absolute value of the criterion for a single model 
has no meaning, and only the relative values between different models are interesting. Considering several candidate 
models, the one that minimizes the AIC (or BIC) is usually considered the best. Comparing with the best one, the 
difference in AIC (or BIC) for other one model is expressed as AAIC = Ax^ij„ + 2Aif (or ABIC = A y^,- -I- AA" In n) . 
Thus one can assess the strength of the models. The rules for judging the AIC model selection are [42 1: when < 
AAICi< 2 model i has almost the same support from the data as the best model, for 2 < AAICi< 4, model i is 
supported considerably less, and with AAICi> 10 model i is practically irrelevant. For BIC, one has: a ABIC of more 
than 2 (or 6) relative to the best one is considered unsupported (or strongly unsupported ) from observational 
data l^ai. 

For ACDM and GCG model in the flat and non-flat universe, the IC values against the model are listed in table 
IIVI According to the table, AIC shows that the flat and non-flat GCG model are supported considerably less by 
current observational data. For BIC selection method, it seems that a more complex GCG model is not necessary to 
explain the current data. In addition, we can see that the current observations much support a flat-geometry universe. 
According to Ref. 4J|, one knows that in the limit of large data points (Inn > 2), AIC tends to favor models with 
more parameters while BIC tends to penalize them. Here it should be noticed that for the used data in our analysis. 
Inn = 6.129, so for BIC the punishment on much-parameter models is more stern. 



^ For small sample sizes, there is a corrected version of the AIC, AICc=AIC+2A'(A' — 1) /(n — X — 1) 41] , which is important for n/ K < 40. 
Obviously, in our case, this correction is negligible. 
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5. Conclusion 

The constraints on the flat and non-flat GCG model as the unification of dark matter and dark energy are studied in 
this paper by using the latest observational data: the Constitution dataset including 397 SNIa, the Hubble parameter 
data, the cluster X-ray gas mass fraction, the baryon acoustic oscillation and the five-year WMAP data. The constraint 
on GCG model parameters are more stringent than the previous papers [l^ j4o|. According to the constraint results, 
since the best fit values of parameters a and Qk are near to zero, it seems that the current observations tends to 
make the GCG model reduce to the flat ACDM model. Furthermore, according to the IC, we can get the same result. 
In addition, we also make a stringent constraint on ACDM model, and in a flat ACDM model it is shown that the 
cosmic age is about, tageiGyr) — 13.7251q'°4'^ (Ict) io.ies (2f )i for using a tophat prior as 10 Gyr < tage < 20 Gyr in 
our calculation. At last, for the CBF constraint method, as a reference we list the best-fit values of the parameters 
in fgas- K = 0.9919, = 0.2089,7 = 1.0299, 6o 0.7728, at = -0.0582, sq = 0.1656, a, = 0.1128 for the flat universe 
and K = 0.9871, rj = 0.2114, 7 = 1.0507, 5o = 0.7749, ab = -0.0950, sq = 0.1741, = 0.0194 in the non-flat case. 
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